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ABSTRACT 



We present sub-milliarcseond resolution imaging and modeling of two nearby 
rapid rotators a Cephei and a Ophiuchi, obtained with the CHARA array - the 
largest optical/IR interferometer in the world. Incorporating a gravity darken- 
ing model, we are able to determine the inclination, the polar and equatorial 
'p_i! radius and temperature, as well as the fractional rotation speed of the two stars 

Q \ with unprecedented precision. The polar and equatorial regions of the two stars 

have ~2000K temperature gradient, causing their apparent temperatures and 
^ \ luminosities to be dependent on their viewing angles. Our modeling allow us to 

determine the true effective temperatures and luminosities of a Cep and a Oph, 
permitting us to investigate their true locations on the H-R diagram. These prop- 
erties in turn give us estimates of the masses and ages of the two stars within a 
^ \ few percent of error using stellar evolution models. Also, based on our gravity 

CN \ darkening modeling, we propose a new method to estimate the masses of single 

^ \ stars in a more direct way through Vsini measurements and precise geometrical 

Q"^ \ constraint. Lastly, we investigate the degeneracy between the inclination and 

^ ' the gravity darkening coefficient, which especially affects the modeling of a Oph. 

Although incorporating Vsini has lifted the degeneracy to some extent, higher 
^ ' resolution observations are still needed to further constrain the parameters inde- 

o3 ' pendently. 
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Introduction 



In the past few years, optical interferometers have resolved the elongated photospheres 
of rapidly-rotating stars for the first time. The emergence of these high angular resolution 
observations of hot stars has shined a spotlight on critical areas of stellar evolution and 
basic astrophysics that demand our attention. For decades, stellar rotation was generally 
overlooked in stellar models and was regarded to have a trivial influence on stell ar evolution 
because most stars are slow rotators, such as the Sun (jMaeder &: Meynetl I2OOOI ) . Although 
the effects of rotation on solar type stars are indeed relatively mild, they are more prominent 
on hot stars. Studies have shown that a l arge fraction of hot s t ars are rapid ro tators with 
rotational velocities more than 120 kms~^( Abt fc Morrelllll995l : lAbt et al.l 120021 ). Virtually 
all the er nission-line B (Be) stars are rapid rotators with rotational velocities of ~ 90% of 
breakup (IFremat et al.ll2005l ). Stars that are rapidly rotating have many unique characteris- 
tics. The centrifugal force from rapid rotation distorts their photospheres and causes them to 
be oblate. This distortion causes their surface brightness and T^ff to vary with latitude, and 
their equatorial temperatures are predicted to b e much cooler than their polar temperatures, 
a phenomenon known as "Gravity Darkening" (jvon Zeipellll924al Jbl). Recent stellar models 
that took r otation into accoun t showed that rapid rotation also affec ts stars' luminosity, 
aburidance ( Pinsonneaultl 1997), evolution, and in creases their lifetime ( Kiziloglu Sz Civelek 



1996; Talon et al 



I997I: I Meynet fc Maederl I2OOOI ) . It is also linked to stellar wind, mass 



MacFadyen et al. 



2001 



loss (e.g.. iMaeder et al.ll2007D. and eyen Ga mma- Ray bursts (IMacFadyen &: Woosleylll999 



Burrows et al.l 120071 ). 



The development of long baseline optical interferometry in recent years has evoked ob- 
servations on several nea rby rapid rotators, fo r instance, Altair, Vega, Achernar, Alderamin 



a Cephei) and Regulus (|van Belle et al 
Domiciano de Souza et al.ll2003 



2006 



2001 



:^Aufde nberg et al. 2006 : Peterson et al. 



van 



Belle et al..i2006i : lMcAhster et al. 



2005 



2006 



Kervella fc Domiciano de Sou: 



Monnier et al.l 120071 ). These studies confirmed the general picture of von Zeipel's grav- 



ity darkening law, but also raised discrepancies between observations and the widely adopted 
standard von Zeipel model (i.e., Te// oc g^jf, where (3 is the gravity darken ing coefficient, 
and = 0.25 for fully radiative envelopes). Particularly, the recent study of iMonnier et al. 



(I2OO7I ) on Altair showed that their model prefers a non-standard gravity darkening law. 
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What is more interesting is that they reconstructed a model-independent image for Altair 
and found a darker-than-expected equator compared to the model. This suggests for the 
first time from observations that the standard gravity darkening law may work only at a 
basic level and other mechanisms need to be introduced to account for the extra darkening. 
To address this issue, we will need more detailed studies and model-independent images of 
rapid rotators. 

In this paper, we present our study of the two nearby rapid rotator a Cephei and 
a. Ophiuchi, observed with the CHARA long baseline optical/IR interferometer array and 
the MIRC beam combiner. The star a Cephei (a Cep , Alderamin, HR 8162, V=2A6, 
H=2.13, d=14.96pc) is the eighth nearest A star in the sky. It was classified as an A7 
IV-V star i n earl y studies, but was recently classified as an A8V main sequence star by 
Gray et al.l (120031 ) . I t is one of the few A stars (including Altair) that are found to have chro- 
mosphere activities (jWalter et al.lll995l : ISimon fc Landsman! 119971 : ISimon et al.ll2002l ). The 



Vsin i measurements of a Cep show large scatter, spanning from ~ 180km s ^ to ~ 245km s ^ 



teernacca fc Perinottolll973l:IUesuri fc Fukudala97ol : iRover et al.ll2007l : lAbt fc Morrd3ll995h . 
Recently, Ivan Belle et all fl200fih studied a Cep using the CHARA array and found it is ro- 
tating close to break-up, and its photosphere is elongated due to rapid rotation. 

The star a Ophiuchi (a Op h , Rasalhague, HR 6556, V = 2.09, if =1.66, d=14.68pc) is a 



nearby subgiant binary system (jWagmanlll946l : iLippincott &: Wagmanlll966l ). and is the sev- 
enth nearest A star in the sky. The primary is a A 5IV sub-giant which was fir st id entified as a 
class H I star but was later corrected to class IV by lAugensen fc Heintz Jl992h andlGrav et al.l 



(I2OOII). Several groups ha v e tried to study th e orbit of the systern (iMcA 



ister fc Hartkopi 



Gatewood 



20051 . 



1984J : iKamper et all Il989l : iMason et al.l Il999l : lAugensen fc Heintzl Il992l : 
etc.), and it was lately determined to have a period of ~ 8.6 yrs and a semi-major axis 
between 0.4" - 0.5". The mass determination of the primary has large scatter, rangin g 



from2M0 to4.9M0 fe.g. lKamper et al.lll989l : lAueensen fc Heintzlll992l : lGatewoodll2005r ). 



The co mpanion, which is approximately a K2y sta r , is thought to have a mass of 0.5- 
1.2 My teamper et al.lll989l : lAugensen &: Heintg|l992l: iGatewoodI 120051). an d is observed to 



be 3.5 mag fainter than the primary in the K band (iBoccaletti et al 



primary was estimated to be ~ 1.6 - 1.7R0( Barnes et al.lll97^ 



rotational velocity Vsini ranges from 210km s ^to 240km s ^ ( iBernacca &: Perinotto 



20011). The size o f the 
Blackwell et al.lll98oh. Its 



1973 



Uesugi fc FukudJll970l : lAbt fc Morrd3ll995l : iRover et al.ll2002h . implying a Oph is spinning 
at a significant fraction of its break-up speed of ~ 270km s^^. 

This paper is organized as follows. We report our observations and data reduction 
schemes in ^ We discuss our aperture synthesis imaging for a Cep and a Oph in ^ and 
present gravity darkening models for both of them in ^ In ^ we present their temperatures. 
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luminosities, and their locations on the H-R diagram. Based on our modehng, we propose a 
new method to estimate the mass of a star in §61 Finally, we discuss our results in ^ and 
present our conclusions in |JH1 



Observations and data reduction 



Our observations were conducted at the Georgia State University (GSU) Center for 
High Angular Resolution Astronomy (CHARA) interferometer array along with the MIRC 
combiner. The CHARA array, located on Mt. Wilson and consisti ng of six 1-meter tele 



scope s, is the longest optical/IR interferometer array in the world (jten Brummelaar et al. 



20051 ). The array is arranged in a Y-shaped configuration and has 15 baselines ranging from 
34m to 331m, providing resolutions up to ~0.5 mas at the H band and ~ 0.7mas at the K 
band. 

The Michigan Infra-Red Combiner (MIRC) was used here to combine 4 CHARA tele- 
scopes together for true interferometric imaging in the H band, providing 6 visibilities, 4 
closure phases and 4 trip le amplitudes simultaneously in 8 narrow spectral channels (see 
Monnier et al.ll2004j . |2006| . for details). MIRC is designed for stable calibrations and precise 
closure phase measurements. It uses single mode fibers to spatially filter the light coming 
from the CHARA beams. The fibers are brought together by a V-groove array in a non- 
redundant pattern. The outgoing fiber beams are then coUimated by a lenslet array and 
are focused by a spherical mirror to form an interference pattern, which consists of six over- 
lapping fringes with non-redundant spatial frequencies. The fringes are focused again by a 
cylindrical lens into a "line" of fringes and are dispersed by low spectral resolution prisms 



with R ~ 50. The dispersed fringes are finally detected by a PIC NIC camera, w. 



onto 8 spectral channels spanning the H band (A=1.5 - 1.8 /xm) (IMonnier et al 



lere they fa ll 



2004 



2006). 



A det ailed description of the control system and software can be found in iPedretti et al. 
J2009h . 



The system visibilities of MIRC are very stable due to our use of single mode fibers. 
However, the atmospheric turbulence changes faster than the 5.5ms readout speed of the 
camera, causing decoherence of the fringes that needs to be calibrated. We therefore observe 
several calibrators adjacent to our targets over each observing night. For the purpose of bias 
subtraction and flux calibration, each set of fringe data is bracketed with measurements of 
background (i.e., data taken with all beams closed), shutter sequences (i.e., data taken with 
only one beam open at a time to estimate the amount of light coming from each beana) , an d 
foreground (i.e., data taken with all beams open but without fringes) (IPedretti et al.ll2009l ). 
Each object is observed for multiple sets. During the period of taking fringe data, a group- 
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delay fringe tracker is used to track the fringes (iThureau et al.ll2006l ). In order to track the 
flux coupled into each beam in "real time" to improve the visibility measurements, we use 
spinning choppers to temporally modulate the light going into each fiber simultaneously with 
fringe measurements. The chopper speeds were set to 25Hz, 30Hz, 35Hz and 40Hz in 2006 
and were increased to 55Hz, 65Hz, 75Hz and 85Hz in 2007 to avoid overlap of modulating 
frequencies caused by chopper drifts. 

We observed a Cep on 4 nights in 2006 and observed a Oph on 8 nights in 2006 and 
2007, using various array configurations optimized for equal Fourier coverage in all directions 
for good imaging. The detailed log of our observations is listed in Table [TJ Figure [1] shows 
the overall baseline coverage of our observations of a Cep and a Oph. 



The data reduction process follows the pipeline outlined by lMonnier et al.l (120071 ). which 
was validated using data on the binary i Peg. In brief, after frame-coadding, background 
subtraction and Fourier transformation of the raw data, fringe amplitudes and phases are 
used to form squared-visibilities and triple products. Raw squared-visibilities are then esti- 
mated from the power spectrum after foreground bias subtraction. After the fiber coupling 
efficiencies are estimated using either the chopping signal or direct fit to the fiber profiles, 
we obtain uncalibrated squared-visibilities and complex triple amplitudes. Finally, calibra- 
tors with known sizes are used to calibrate the drifts in overall system response before we 
obtain the calibrated squared-visibilites, closure phases, and complex triple amplitudes. The 
adopted sizes of our calibrators are listed in Table [2J Corresponding errors of the data are 
estimated by combining both the scatter of the data and calibration errors. 



3. Aperture Synthesis Imaging 



We employed the pu blicly-available appl ication "Markov-Chain Imager for Optical In- 
terferometry (MACIM)" (jlreland et al.ll2006l ) to reconstruct image s for a Cep and a Oph. 
The application applies the Maximum Entropy Method (MEM) (INarayan fc Nityananda 
1986[) widelv used i n radio synthesis imaging, and has been vahdated on other test data 
( jhawson et al. 2006). Since the photosphere of a star has a sharp emission cut-off at the 
edge, which is imprinted in the highest spatial frequencies that cannot be observed, we con- 
strain the field of view of the images within an ellipse to avoid spreading-out of the fiux by 
the MEM procedure at the edge of the star. This constraint is appropriate for a Cep and 
a Oph due to their lack of any circumst ellar emission outside of their photospheres. The 
details of this approach can be found in iMonnier et al.l (120071 ). The ellipse prior is found 
by conducting MACIM imaging on a grid of ~ 400 different ellipses with uniform surface 
brightness, spanning a range of possible sizes, axial ratios, and position angles. To ensure 
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the smoothness of the image, we also de-weighted the high resolution data with a gaussian 
beam of 0.3 milliarcsec FHWM, an approach usually applied in radio synthesis imaging. 
The image with the global maximum entropy is then taken as the final result. We treated 
each wavelength channel as providing a distinct set of (u, v) plane coverage, ignoring any 
wavelength-dependence of the image itself. This assumption is well justified for a Cep and 
a Oph since the brightness profiles of their photospheres are almost identical in all channels 
in the H band. 

Figure [2] shows the reconstructed image of a Cep [xv' = 1-10). Its photosphere is 
well resolved and appears elongated along the east-west direction. The bright region at the 
bottom with Teff above 7000K (left panel) is later identified close to the pole and the dark 
belt below 6500K is the equator - a direct confirmation of the gravity darkening effect. The 
image implies the pole of a Cep is medium inclined. The very top of the image becomes 
bright again since the photosphere is brighter toward the poles. The right panel of Fig|2] 
shows the orientation of a Cep based on the model in [|H It shows that the bright spot in 
the image is in fact above the pole as the pole of a Cep is limb-darkened. The squared- 
visibilities, closure phases, and triple amplitudes derived from the image are compared with 
the data in Figure [31 SJ and O 

Although we have tried intensively to reconstruct an image for a Oph, we are unable to 
find a reliable solution for it. This is because the brightness distribution of a stellar surface 
is mainly imprinted in our closure phases. The closure phase is only sensitive to asymmetric 
structures of the object, while a symmetric object only gives either 0° or 180° closure phases. 
The squared-visibilities of our data are less constraining due to their relatively large errors. 
The near equator-on inclination of a Oph (see §4.2p makes its brightness distribution nearly 
symmetric, providing too few non-zero closure phase signatures to constrain the image. 
Therefore, we could not obtain a reliable solution for a Op h in the image reconstruction. 



We have also pursued other imaging programs such as MIRA (jThiebautll2008l ). and obtained 
similar results in our preliminary efforts (Thiebaut 2008, private communication). Thus we 
only present the model of a Oph in this paper. As we will see in §4.2[ the lack of non-zero 
closure phase signatures of a Oph also brings similar issues to our modeling, causing high 
degeneracy to the inclination and the gravity darkening coefficient. 



4. Surface Brightness Modeling 



In addition to synthesis imaging, we construc t rapid rotator models to fit the data of 
both stars, following the prescription described in lAufdenberg et al.l (120061 ) and references 
therein. Specifically, we assume a Roche potential (point mass) and solid body rotation in 
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our model, and use the von Zeipel gravity darkening law ( von Zeipel 1924a| jbl) to characterize 
the latitudinal temperature profile. Six parameters are used to define the models, including 
the stellar radius and temperature at the pole, the angular rotation rate as a fraction of 
breakup (w), the gravity darkening coefficient the inclination angle, and the position 
angle (east of north) of the star. To ensure accuracy of the models, we construct them at 
four different wavelength channels across the H band. The intensity and limb darkening at 
each point of t he stellar surface is interpolated using the stellar atmosphere models of Kurucz 
( lKuruczlll993l ) as a function of local temperature, gravity, viewing angle, and wavelength. 
The 3D surfaces of the models are generated using patches with uniform surface areas to 
avoid over-sampling at the poles or under-sampling at the equators, and also to speed up 
the computation. A direct Fourier transform is then used to convert the projected intensity 
model to squared- visibilities, closure phases and triple amplitudes 0. In addition, we also 
force our model to match the V and H band photometric fluxes obtained from the literature 
(see Tables |3l Hj) to constrain the temperature range. 



4.1. a Cep 



We first fit the data of a Cep with the standard von Zeipel gravity darkening model 
for fully radiative envelopes (i.e., Tg// oc S'fjj, where {3 = 0.25; hereafter, the standard 
model). The Levenberg-Marquardt algorithm is applied for the least-square minimiza- 
tion and the pararneter s paces are extensively se arched in the fit . We assume M = 2.0 
MfT^ dvan Belle et al.l 120061). dista nce = 14.96 pc (jPerryman et al.l Il997l ). and metallicity 
[Fe/H] = 0.09 (IGray et al.ll2003l ) in the model. The left panel of Figure [6] shows the best-fit 
standard model of a Cep, with an overall goodness of fit xl of 1.21. The model shows the 
photosphere of a Cep is elongated, with a bright polar region at the bottom and a dark 
equator above it - generally consistent with the synthesized image in Fig|2l Our standard 
model yields an inclination of 64?9 ± 4? 1 and a position angle of — 178?3±4?1, consistent with 
the ellipse fit of Ivan Belle et al.l (l2006l . hereafter VB06), which gave a position angle of -177° 
(or 3° depending on the definition). However, both the inclination and the position angle 
of their gravity darkening model {i = 88?2, P.A. = 17° or — 163") differ from our results, 
as we have better UV coverage and also closure phase information which is very sensitive 



^We have validated our model by comparing with anot her independent model from Jason Auf- 
denberg (private communication) on the data of Vega from Aufdenberg et al. ( 2006h . We also com- 
pared the model using Kurucz limb darkening with one using PHOENIX limb darkening and found 
the difference is negligible. The data and models we used for the comparison are available at 
,http:/ /www.astro. lsa.umich.edu/ ^mingzhao/rapidrot.phfi 
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to asymmetric structures. Our model indicates a Cep is rotating very fast, at 92.6% of 
its break-up speed. The temperature at the poles is ~2400K higher than at the equator, 
while its radius at the equator is 26% larger than at the poles. The best-fit parameters of 
the standard model are listed in the second column of Table [31 Since the calibration errors 
vary from night to night, we estimate the parameter errors by bootstrapping the data from 
different nights (i.e., treat each night of data as a whole and randomly sample all of the 
nights with replacement, so that the correlations of data within each night can be taken into 
account ) and fitting the parameters to the resampled data. We then iterated this procedure 
hundreds of times. 

In addition to our data, we also combine the squared-visibilities from VB06 (here after 
"Classic data") into our fit. The combined fit gives a slightly higher inclination, but all 
parameters are still consistent with our original fit. The total Xv^ of the combined fit is 
1.25. However, the Xv^ of the Classic data (x;^^=2.0) is very large although it is slightly 
better than the original result of VB06 {xJ^ =2.16), implying that either the Classic data 
have additional un-calibrated errors or the model needs more degrees of freedom. We first 
look into a free /3 in the model. Indeed, the von Zeipel theory suggests that the standard 
gravity darkening coefficient (/? = 0.25) only applies to pure radiative envel opes. However , 



it is uncertain if a Cep is pure radiative or not. The atmosphere models of iKuruca (Il979l ) 
suggest that, for an atmosphere with Tg// > 75007^^ and log (7 ~ 4, like the polar areas of 
a Cep, convection should have very little or no effect. But it starts to play a role when 
temperature and log q drop below those numbers. In addition, the evolution models of (5 



calculated by IClaretl (Il998l . 120001 ) also indicate that, for a 2-2.5 Mq star, convection starts 
to take place once Te// is below ~7900K. For the case of a Cep, although its Te//s at the 
polar higher than 8000K, they drop to only ~6700K in the equator, implying that 

convection may have effects in the equatorial areas and /3 may deviate from the standard 
value. Therefore as a preliminary effort, we extend the standard von Zeipel law to a free /3. 

The new combined /5-free fit gives a Xv^ of 2.11 to the Classic data, similar to the original 
VB06's result. But it prefers a /5 of 0.22 rather than the 0.08 value of VB06. To address 
this issue, we tried to fit the combined data at a fixed /5 of 0.08 instead, but only obtained a 
total Xv^ of ~ 6.5, much worse than the previous result. In addition, we also fit the Classic 
data only but found it is too hard to constrain the model due to the small amount of data 
and lack of phase information. Therefore, due to possible uncertainties of the Classic data, 
we applied the /?-free model to the MIRC measurements only, and the results are shown in 
the third column of Table [31 The best-fit model is shown in the right panel of Figure [51 The 
squared-visibilities, closure phases and triple amplitudes of the /3-free model are compared 
with the data in Figures [31 [H and [5l respectively. 
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The right panel of Figure [6] shows that the /3-free model is more consistent with the 
synthesized image in Figj2]than the standard model. The xJ^ of closure phase is significantly 
improved in the new best-fit although the xJ^ of the triple amplitude is slightly larger. 
Figure [7] illustrates the space of inclination and /? for a Cep, showing the value of (3 is 
well constrained in the new model and is slightly lower than the standard value of 0.25. We 
also test the corresponding Vsini of the models in FigjTl The peak of the x^ space falls 
inside the green box, consistent with the observed range of Vsini. The new model prefers a 
lower inclination of 55?70 ± 6?23, a higher rotational speed of 94% of break-up, and a similar 
position angle. The new best-fit temperatures at the poles and the equator are both cooler 
than those of the previous standard model. 

In addition to using an average (3 throughout the stellar surface as applied above, we 
are also pursuing fitting /3 as a function of latitude. This approach will be presented in a 
future work with higher resolution data. 



4.2. a Oph 



We also start with the standard gravity darkening m odel (/3=0.25) fo r a Oph. We 
assume mass = 2.10 M© (see §5]) and dis tance = 14.68 pc ( Gatewoodl l2005l ) in the model. 
The metallicity [Fe/H] of a Oph is -0.16 (lErspamer &: Northll2003l ). thus a Kurucz grid with 
metallicity of -0.2 is applied. Figure [8] shows the best-fit standard model of a Oph. The 
best-fit parameters are listed in Table |H The associated errors of the parameters are also 
obtained using the bootstrap procedure described in §4.11 The squared-visibihties, closure 
phases and triple amplitudes of the model are compared with the data in Figures [HI [IHl and 
[TTl respectively. The model shows that the photosphere of a Oph is also elongated and has 
two bright polar areas and a dark equator. Its radius at the equator is ~ 20% larger than 
at the poles. It is seen nearly equator-on with an inclination of 87?70 ± 0?43. The model 
also shows that a Oph is rotating at 88.5% of its break-up speed and the poles are ~ 1840K 
hotter than the equator. 

In the standard model, the Xv"^ of the closure phase only reaches 1.33 (Table Hj), sug- 
gesting that we may need extra degrees of freedom to improve the fit. Therefore, following 
our approach for a Cep, we extend the standard model of a Oph to a free /?. However, 
although we have searched the parameter space extensively, we cannot find a unique /3-free 
model for a Oph due to the same reason that we encountered in imaging. As we mentioned 
in ^ this issue stems from the near equator-on and symmetric brightness distribution of a 
Oph, causing the closure phases to be mostly 0° or ±180° (as shown in Fig JTOj) and hence 
lack of enough non-zero signatures to constrain the model when (3 is free. 
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Figure [12] shows the xl space of inchnation and (3 for a Oph. Unhke the single peak of 
a Cep, a Oph has several peaks spreading over a large range of inclination and (3, indicating 
the inclination and jS are highly degenerate and suggesting it is difficult to constrain a 
unique /5-free model. Nevertheless, the corresponding Vsini values around the largest peak 
at /3 ~ 0.08 fall outside the observed range of 210 - 240 km s~^ (enclosed by the green box 
in Figure [12]), suggesting the peak is not real but only due to the degeneracy of /3 and 
inclination. In addition, the peak around (3 ~ 0.08 corresponds to a fully convective star 
according to iLucy] (119671 ). But it is unlikely for an A5 star to be fully convective, especially 
when its polar temperature is as high as 9300K. Therefore, we can rule out t he large s t peak 
around (3 ~ 0.08. Furthermore, the gravity darkening evolution models of I Claret! (120001 ) 
show that the value of (3 should be much larger than 0.15 for a ~ 2M0 star with average 
Teff higher than 7500K, like a Op h. The s econd peak around /5 ~ 0.15 in FigJT^ however, 
is not consistent with the models of IClaretl (120001 ) although is inside the Vsini range. Thus, 
in this study we still prefer the other peak around the standard /? = 0.25 model for a Oph. 
To break down the degeneracy and constrain the value of f3 more accurately, we will need 
more observations with higher resolution, especially in the visible where limb- darkening and 
gravity darkening are more prominent. 



5. Physical properties and comparison with stellar evolution tracks 



In addition to the model parameters, we also calculate the true and apparent effective 
temperatures and luminosities for the two stars in Table [3] & [U The true luminosity is 
estimated by integrating local aTf,ff{9Y (where a is the Stefan-Boltsman constant) over 
the stellar surface, and the true Tg// is estimated from the total luminosity and the total 
surface area of the star. The apparent luminosity is obtained from L = 4:7id'^Ff,oi, where 
the bolometric flux Ftoi is calculated by integrating the specific intensity over the whole 
spectrum and the projected angular area of the star. The apparent temperature is obtained 
nd^Fhoi/Af 



from crT^jj- 



'■proj, where Aproj is the projected area. 



The true Teff and luminosity of a Cep are very close to its apparent values due to its 
medium inclination (see Table [3]). Its true Teff from the /5-free model is 7510 =fc 160K, c lose 
to although slightly cooler t han the ~7700K estirnate of VB06 and [Gray et al.l (120031 ) . as 
well as the 7740K estimate of Malagnini fc Morossi ( 199oh. Its true l umin osity is 18.1 ± 1.8 
Lq, consist ent with the 17 L^t^ estima te from iMalagnini fc Morossi] (ll990l ) and the 17.3 Lq 
estimate of ISimon &: Landsman] (11997] ) 



The deviation of a Oph's true Tg// and luminosity from its apparent values is very signif- 
icant because of its near equator-on inclination. Its true Teff from the standard model is esti- 
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mated to be 8250±100K. Its apparent Tg//, on the other hand, is 7 950K based on the model, 
consis tent with the apparent value of 788 3 =fc 63 K calculated b y iBlackwell &: Lynas-Gray 
f|l998l ) and the value of 8030 ± 160 K by iMalagnini fc Morossil Jl99oh. Its app arent lumi- 
nosity is 24.3 Lq, in agreement with the 25.1 Lq value of IMalagnini &: Morossil (Il990l ) but 
smaller than its true luminosity of 30.2 ± 1.3 Lq. 

Because rapid rotators are hotter at the poles and cooler at the equators, their apparent 
temperatures are therefore dependent on their inclinations, which can easily introduce large 
biases to the observed values. To investigate this effect, we plot in Figure [T5] the differences 
between the true and apparent values of Te//s and luminosities as a function of inclination, 
scaled with their true values. The plots show that when a star is inclined by ~ 54°, its 
apparent Te// and luminosity seen by the ob servers will b e equa l to their true values, just as 



the case of a Cep and similar to the result of iGillich et a 



on, such as Vega (lAufdenberg et al.ll2006l : iPeterson et al. 
exceed the true value by 



(120081) . When the star is seen pole 
2006h . its apparent temperature can 
5%, and the luminosity can exceed by ~ 40 — 50% or even larger 



depending on the speed of the rotation, which explains the reaso n that Vega's lurn i nosity 
was largely overest i mated for a long time until recent studies of lAufdenberg et al.l (l2006l ) 
and iPeterson et al.l (120061 ). On the other hand, when a rapid rotator is equator-on, as the 
case of a Oph, its apparent temperature and luminosity can be underestimated by ~ 4% 
and ~ 20% respectively. The rotation speed of the star also affects the differences between 
its true and apparent values - the faster the star rotates, the larger the difference we see. 

Our estimates of the true Te//S and luminosities of a Cep and a Oph also allow us to 
understand their current evolutionary status better. In Figure [14] w e plot the H-R diagrarn 
and the corresponding Y'^ stellar evolution tracks and isochrones (jPemarque et al.l |200J) 
for a Cep and a O ph. Their possible ranges of locations on the H-R diagram (also called 
"inclination curve" , iGillich et al.ll2008l ) are also shown in the plots. The top panel shows that 



a Cep appears to be an A9 type star on the H-R diagram based on its ap parent temperatur e 
and luminosity (filled triangle). However, it is classified as an A8V star by lGray et al.l (120031 ) . 
earlier than that inferred from the top panel. Similarly, in the bottom panel o f Figure [T^ 



a Oph appears roughly as an A6.5 type star. Its apparent spectral type from iGray et al. 



(I2OOII ) is A5IV, also earlier than that inferred from the figure. We infer that this is because 
the spectra of the two stars are dominated by spectral lines from the hotter and brighter 
polar regions, causing their overall spectral classification to be biased toward the types of 
their poles which appear earlier than other regions of the stars. Therefore, for the case of an 
equator-on star, such as a Oph, although its apparent effective temperature is lower than 
its true temperature due to the inclination, its spectral type derived from spectroscopy can 
compensate this effect and make it look closer to its true spectral type. However, for a 
pole-on star such as Vega, this bias can not be compensated, and the spectral types derived 
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from both spectroscopy and apparent temperature will appear earlier than its true type. 
This phenomenon indicates that the spectral types of rapid rotators are not only biased by 
their inclinations, but also by the spectral lines of their polar regions. 

Using the Y"^ models, we estimate that a Cep has a mass of 1.92 ± 0.04 Mq, slightly 
smaller than the estimate of VB06. Its age is estimated to be 0.99 ± 0.07 Gyrs. We also 
estimate that a Oph has a mass of 2.10±0.02 Mq, and an age of 0.77±0.03 Gyrs. Its apparent 



position in the H-R diagram, howeve r, indicates a lower rnass of 1.99 
consis tent with the 2.0 estimate of iMalagnini fc Morossil (119901 ) and 



Vlf;^, which is again 



Augensen fc Heintz 



(I1992I ). Howev er, this value is much lower than the 2.84 M0 value of iGatewood J2005h and 
the 4.9 M0 of iKamper et al.l (119891 ). To address the differences, we derive the mass range 
of a Oph usi ng our new meth od of estimating ma s s in t he next section and conclude 
the result of iGatewoodI (120051 ) and iKamper et al.l (119891 ) can be ruled out. The estimated 
masses and ages of a Cep and a Oph are included in Tables [3] and H] respectively. 

We note that the Y"^ models are for non-rotating stars, whereas both a Cep and 
a Oph are rapid rotators . The fact that rotation may extend the main-sequence lifetime 
(IKiziloglu fc Civeleklll996l : iMaeder fc Meynetll2000l ) implies that our age estimates may not 
be accurate and needs further investigation. We also note that the masses of a Cep and 
a Oph are both estimated based on non-a-enhanced models. Studies ha ve shown that 
rapid rotation can c hange the abund ance of a star (e.g., iPinsonneaultl 119971 ) and enhance 



the a-rich elements ( lYoon et al.ll2008l ). resulting in very different estimates of its mass and 
age. Hence to derive the masses of a Cep and a Oph more accurately, detailed abundance 
studies are required to determine if they are a-enhanced and what abundance to use for 
their evolutionary models. 



A new method to estimate the mass of a star 



Mass is the most fundamental property of stars. The deter mination of stella r masses 

2007h. stellar 



mostly relies on orbital measurements of binary systems (e.g., IZhao et al 



evolu tion models together with measurements of other stellar properties fe.g.. lyan Belle et al. 



2006 ). and asteroseismology together with measurements of stellar radii (e.g.. ICreevey et al. 



20071 ). Here we propose a new method to estimate the mass of a star based on our modeling 
of rapid rotators. 

Since we can determine the inclination, equatorial radius and the fractional rotation 
speed of a rapid rotator from our model, we therefore can combine the model of a rapid rota- 
tor with its mass to estimate the equatorial velocity and the Vsin i value. We can also reverse 
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the process, taking a precise measurement of Vsini and a best-fit rotator model to determine 
the mass of a star. This approach is most suitable for radiative rapid rotators which can be 
interpreted by the standard gravity darkening model, and also non-fully-radiative rotators if 



a more sophisticated fiuid model is co nstructed (e.g.. IJackson et al.ll2004j : iMacGregor et al. 



20071 : lEspinosa Lara fc Rieutordll2007l ). For stars with less accurate models, we can also use 



this method to roughly estimate their masses. The precision of Vsini is also crucial for a 
precise mass estimate. As a preliminary test, we first apply this method to a Cep and a 
Oph. 

The Vsini range of a Cep (180 kms~^- 245 kms~^) corresponds to a large mass range 
of 1.3 to 2.4 M0 based on the /?-free model in §4.11 The mass of a Cep determined from 
stellar models, on the other hand, is 1.92±0.04 Mq (see well within the mass range given 
by Vsini. Similarly, the Vsini range of a Oph (210 kms~^- 240 kms~^) gives a mass range of 
1.7 M0 to 2.2 M0 when combined with the model in §4.2[ Its mass determined f rom stellar 
model s, 2.1 ±0.02 Mfn (se e JjSl), is also within the range. By contrast, the study of lGatewood 
( I2OO5I ) and lKamper et al.l (119891 ) gave a mass of 2.84 M0 and 4.9 M0 to a Oph respectively, 
far outside the range given by Vsini, and hence can be ruled out. Since a Oph is also a 
known astrometric binary, it is the ideal target to further test this new method by comparing 
its mass with that determined from the astrometric orbit. We are currently pursuing this 
study (Oppenheimer et al.2008, private communication) and will also present it in a future 
work. 



Discussion 



Although the /3-free model of a Cep is consistent with the synthesized image (FigI2]) in 
basic features such as the bright pole and the dark equator, we also notice that the equator 
of the image is darker and cooler tha n that of the model - a phenomenon seen in a previous 
study of Altair (jMonnier et al. 2007). The existence of the darker-than-expected equator on 
both stars implies that the extra gravity darkening may be real. However, it can also be 
due to a systematic effect of the imaging program. To confirm this conclusion we will need 
further studies such as model-independent latitudinal temperature profiles. 

Our models show that both a Cep and a Oph have polar temperatures well above 8000K 
and eq uatorial t empe ratures below 7500K, which means, according to the stellar atmospheric 
grid of iKuruczl (119791 ) . the polar areas of a Cep and a Oph are radiative and their equators 
can have convections, especially for a Cep as its equatorial temperature is lower than that 
of a Oph. Since the existence of convecti on tends to l ower the value of the average gravity 
darkening coefficient f5 of the whole star (jClaretlll998l ). it may be the cause of /? < 0.25 in 
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the P -free model of a Cep. The unusually strong c hromosphere activity of a Cep among A 



stars (jWalter et al.lll995 



Simon fc Landsman! 1 19971 ) also provides evidence to the convective 



layers since the chromosphere is directly linked to magneto-convection. Another A star with 
strong chromosphere activities, Altair, is also a rapid r otator spinning at 92 % of its break- 
up speed and has an equatorial temperature of 6860K (IMonnier et al.l 120071 ) . This suggests 
that although A stars are genera lly considered to h ave no chromospheres due to their very 
thin or lack of convective layers (jSimon et al.ll2002l ). rapid rotators may have exceptions at 
their equators due to gravit y darkening. This is also consiste nt with the conclusion from 
the hydrodynamic model of lEspinosa Lara fc RieutordI (120071 ). This effect may also shed 
some light on the searches for the onset of chromosp here and the tran sition from radiative 
to convective envelopes among early type stars (e.g., 



Simon et al.ll2002h . 



Since convection also tends to smear out the temperature differences between the hot and 
cool regions of the stellar surface a nd make their intensity contrast l ower, other mechanisms 



such as differential rotation (e.g., lEspinosa Lara &: RieutordI 120071 ) may also exist in the 



equators of these stars in order to make the equator darker and cooler as in the image. For 
instance, a faster differentially spinning equator will have stronger gravity darkening, thus 
will appear darker than that of the standard model. However, the darker equator, if it is real, 
can also be caused by a very different form of gravity darkening law. To further address this 
issue, we will need detailed line profile studies and images at visible since gravity darkening 
is more prominent in the visible than in the H band. 

The 87?70 inclination of a Oph differs from its orbital inclinat ion by about 27° {i ~ 
115°, iKamper et al.lll989l : lAugensen fc Heintg|l992l : iGatewoodI 120051 ). indicating the spin of 
a Oph is not coplanar with its orbit. Even more interesting, the orbit of the binary is highly 
eccentric (e ~ 0.8, Kamper et al. 1989 and Gatewood 2005; e = 0.57, Augensen & Heintz 
1992), implying the non-coplanarity and the high eccentricity of the system may be related 
to each other through interactions of the two stars with their disks in their early formation 
stages. 



8. Conclusion 

We have modeled the surface brightness distributions of a Cep and a Oph using the 
gravity darkening model. We have also reconstructed an aperture synthesis image for a 
Cep, but no reliable image for a Oph is available due to its lack of closure phase signatures 
caused by its nearly symmetric brightness distribution. The image of a Cep shows the star 
is oblate and its equator is darker than its poles, directly confirming the gravity-darkening 
phenomenon. The models show that both stars are rotating close to their break-up speed. 
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They both appear oblate and have large latitudinal temperature gradient due to gravity 
darkening. A standard gravity darkening model of (3=0.25 is adopted for a Oph, and its 
inclination is determined to be 87.70°. For a Cep, a (3 = 0.216 model fits the data better 
and also agrees better with the image. It has a medium inclination angle of 55.70°. 

Our models also allow us to calculate and compare the true Te//s and luminosities of 

the two stars with their apparent values. We show that a Oph has a tnie Teff of 8250K 
and luminosity of 30.2 Lq, significantly larger than its apparent values due to its equator-on 
inclination. The true Teff and luminosity of a Cep, on the other hand, appear very close 
to its apparent values because of its medium inclination. The spectral classification of the 
two stars from literatures, however, suggests earlier spectral types for both stars than that 
derived from their apparent Te//S and luminosities. We infer that this is because the spectra 
of the two stars are dominated by lines from their hotter and brighter polar regions which 
appear much earlier in spectral type than the other regions of the stars, causing their overall 
spectral classification to be biased toward their polar areas. 

The temperatures and luminosities in turn allow us to make rough estimates of the 
masses of the two stars through stellar evolution models. The mass of a Cep is estimated to 
be 1.92 M0, and the mass of a Oph is 2.10 Mq. However, due to possible abundance anomaly 
caused by rapid rotation, the exact masses of the two stars still have to be scrutinized when 
a detailed abundance analysis is available. 

Our gravity darkening models also allow us to propose a new method to estimate the 
masses of rapid rotators together with precise measurements of Vsmi. We have tested this 
method on both stars and found our mass estimate from the stellar models are within the 
range. The star a Oph will be a good target to further test this method as it is also an 
astrometric binary. 

Our models show that the equatorial temperatures of a Oph and especially a Cep are 
low enough to meet the onset conditions of convection, implying that convections in the 
equatorial region can be a reason of the unusually high chromosphere activities of a Cep. 
Although the a Cep model agrees with its image in general, the image shows extra darkening 
at the equator which is not expected by our gravity darkening model but is consistent with 
the previous result of Altair. This effect, if is real, is most likely caused by differential 
rotation of the star. But to further confirm the conclusion, detailed high resolution line 
profile analysis and images at visible are needed. 
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Fig. 1. — Baseline coverages for a Cep and a Oph. The longest baselines in the observations 
are 251m and 329m for a Cep and a Oph, corresponding to resolutions of 0.68mas and 
0.52mas respectively. The UV coverage can be obtained by dividing these two plots by 
corresponding wavelengths. 
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Table 1. Observation logs for a Opli and a Cep 



Target 


Obs. Date 


Telescopes 


Calibrators 


Chopper 


ct Oph 


UT 2006Jun20 


W1-W2-S1-S2 


a Sge 


no 




UT 2006Jun21 


W1-W2-S1-S2 


C Oph, 7 Ser 


no 




UT 2006Aug28 


S2-E2-W1-W2 


11 Peg 


no 




UT 2006Aug29 


S2-E2-W1-W2 


7 Lyr, v Peg 


no 




UT 2006Aug30 


S2-E2-W1-W2 


7 Lyr 


yes 




UT 2006Aug31 


S2-E2-W1-W2 


7 Lyr, v Peg 


yes 




UT 2007MaylO 


S1-E1-W1-W2 


C Oph, T Aql 


yes 




UT 2007Mayl2 


S1-E1-W1-W2 


C Oph, r Aql 


yes 


Q Cep 


UT 2006Oct09 


S2-E2-W1-W2 


29 Peg, V And, C, Per 


yes 




UT 2006Octll 


S2-E2-W1-W2 


V And, C Per 


yes 




UT 2006Octl2 


S2-E2-W1-W2 


29 Peg, C Per 


yes 




UT 2006Octl6 


S2-E2-W1-W2 


29 Peg, V And 


yes 



Table 2. Calibrator diameters 



Calibrator 


UD diameter (mas) 


Reference 


a Sge 


1.32 ±0.02 


Uniform-disk fit to PTI archive data'' 


COph 


0.51 ± 0.05 


Hanburv Brown et al. fl974) 


7 Ser 


1.21± 0.05 


Uniform-disk fit to PTI archive data 


7 Lyr 


0.74 ±0.10 


Loggctt et al. (1986) 


V Peg 


1.01 ± 0.04 


Blackwcll & Lvnas-Grav fl994) 


r Aql 


1.10 ± 0.01 


Merand et al. f2005. 20061 


29 Peg 


1.0 ± 0.1 


MIRC measurement 


V And 


1.17 ± 0.02 


Boden 2008*" 


C Per 


0.67 ± 0.03 


getCaF 



''available at |http : / / mscweb . ipac .caltech.edu / mscdat-pti 
''SED fit, private communication 
ht tp : / / mscweb . ipac . caltech . edu/ gcWeb /gc Web . j sp 
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Table 3. Best-fit and physical parameters of a Cep 



Model Parameters 


Standard (/3 = 0.25) 


Non-standard (/3-free)* 


Inclination (degs) 


64.91 ± 4.11 


55.70 ± 6.23 


Position Anglo (dogs) 


-178.26± 4.10 


-178.84 ± 4.28 


Tpoi (K) 


8863± 260 


8588 ± 300 


Rpoi (R©) 


2.199 ± 0.035 


2.162 ± 0.036 


Te, (K) 


6707 ± 200 


6574 ± 200 


Reij ( Rq) 


2.739± 0.040 


2.740 ± 0.044 




0.926 ± 0.018 


0.941 ±0.020 




0.25 (fixed) 


0.216± 0.021 


Model V Magnitude'' 


2.45 


2.45 


Model H Magnitude'' 


1.92 


1.91 


Model V sini (km/s) 


237 


225 


Total xi 


1.21 


1.18 


Vis2 xl 


0.79 


0.80 


CP xl 


1.43 


1.27 


T3amp xl 


1.71 


1.76 


Other Physical Parameters 


Trne T:^f[(K) 


7690 ± 150 


7510 ± 160 


True Luminosity (Lq) 


20.1 ± 1.6 


18.1 ± 1.8 


Apparent 'T^ff(K) 




7510 


Apparent Luminosity ( Lq ) 




17.9 


Mass (Mq)'= 




1.92 ± 0.04 


Age (Gyrs)'^ 




0.99 ± 0.07 


[Fe/H]'^ 




0.09 


Distance (pc)^ 




14.96 


*The /3-free model is adopted 


as the final model, see text of i|4.1lfor detail. 


''V magnitude from literature: 2.456 ± 0.002 JPerrvman et alJll997^ 



''H magnitude from literature: 2.13 ± 0.18 llCutri et alj|2003h 
'^Based on the stellar evolution model jPemargue et al]|2004l) . 
' iGrav et all ll2003h 
i Ferryman et a~ I lll997h 
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Table 4. Best-fit and physical parameters of a Oph 



Model Parameters 


Standard (/3 = 0.25) 


Inclination (degs) 


87.70 ± 0.43 


Position Angle (degs) 


-53.88 ± 1.23 


Tpo, (K) 


9300 ± 150 


Rpoi (R-o) 


2.390 ± 0.014 


Te, (K) 


7460 ± 100 


Reg (R.q) 


2.871 ± 0.020 




0.885 ± 0.011 


P 


0.25 (fixed) 


Model V Magnitude^" 


2.086 


Model H Magnitude'' 


1.66 


Model V sini (km/s) 


237 


Total xl 


0.91 


CP xl 


1.33 


Vis2 xl 


0.72 


T3amp 


0.81 


Other Physical Parameters 


True T^ff{K) 


8250 ± 100 


True Luminosity ( Lq ) 


30.2 ± 1.3 


Apparent Te^^(_R') 


7950 


Apparent Luminosity (Lq) 


24.3 


Mass {Mq)'= 


2.10 ± 0.02 


Age (Gyrs)'^ 


0.77 ± 0.03 


[Fe/H]<^ 


-0.16 


Distance (pc)"^ 


14.68 



''■V magnitude fro m literature: 2.086 ± 0.003 
jPerrvman et al.|[l997l ) 

''H magnitude from l iterature: 1.66 ili . 03 (weighted 
average of fluxes from: lAlonso et al.|[l998l : ICohen et al.l 
Il999l:lcutri et al.ll2003h 

Based on the stellar evolution model 

jPemargue et al.|[20o3) . 

' lErspamer fc NorthI ||2003|) 

iGatewoodI ||2005|) 
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Fig. 2. — Reconstructed MACIM image of a Cep. The left panel shows the contours of 
local brightness temperature. To help visualize the geometry of a Cep, the right panel shows 
its latitude and longitude using the positions from the standard model discussed in Jgl The 
white circle at the bottom-left corner of the left panel shows the size of the convolving beam 
that we use for the image reconstruction. The total of the image is 1.10. The resolution 
of the image is 0.68 milliarcsec. 
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Fig. 3. — a Cep squared- visibilities from the MACIM image (solid lines) and the gravity 
darkening model (/9 = 0.216, dashed lines) vs. data (filled points with error bars). All four 
nights (2006 Oct09, 11, 12, 16) are shown here. The of the image's squared-visibilities 
is 0.87, while that of the model is 0.80. Each row stands for a different baseline, while the 
columns indicate different times of observation. The eight data points in each panel indicate 
the eight spectral channels of MIRC across the H band. (Please refer to the electronic edition 
if the type size is too small.) 




Fig. 4. — Similar to FigjS] but showing the closure phases for a Cep. The solid lines stand 
for the closure phases of the MACIM image, and the dashed lines stand for the model. Each 
row stands for a different telescope triangle. The of the image's closure phases is 0.95, 
while that of the model is 1.27. 
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Fig. 5. — Similar to FigO] but showing the triple amplitudes for a Cep. The of the 
image's triple amplitude is 1.63, while that of the model is 1.76. 
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Fig. 6. — The gravity darkening models of a Cep . The contours indicate the local brightness 
temperatures on the surface of the star. The left panel shows the best-fit standard gravity 
darkening model (/5 = 0.25) overplotted with the temperature contours from Figure [2l The 
total of the standard model is 1.21. The right panel shows the best-fit /5-free model, also 
overplotted with the temperature contours from Figure [2], and has a total of 1.18. The 
resolution of the data is 0.68 milliarcsec. 
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Fig. 7. — The xl surface of f5 and inclination for a Cep. The corresponding probabihty peaks 
aX (3 ^ 0.22 and i ~ 56°. The black contours show the 1-a, 2-a, and 3-a levels of confidence 
interval, scaled to match the errors of /3 and inclination estimated from bootstraping. The 
area inside the green box indicates the region where the corresponding Vsini values are 
within the observed range of 180 - 245 kms~^. 
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Fig. 8. — The best-fit standard gravity darkening model of a Opli. Tlie contours in tfie left 
panel indicate the local brightness temperatures on the surface of the star. The right panel 
shows the latitude and longitude of a Oph to help visualize its geometry. The resolution of 
the data is 0.52 milliarcsec. The total xl model is 0.91. 
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Fig. 9. — OL Oph squared-visibility model (standard /3 = 0.25, solid lines) vs. data 
(filled points with error bars). Four out of eight nights (2006Jun21, 2006Aug29, 31, and 
2007Mayl2) are shown here. Each row stands for a different baseline, while the columns in- 
dicate different times of observation. The eight data points in each panel indicate the eight 
spectral channels of MIRC across the H band. The total is 0.72 for the squared-visibility 
only. (Please refer to the electronic edition if the type size is too small.) 
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Fig. 10. — Similar to Figl9]but showing the closure phase for a Oph. Each row stands for 
a different telescope triangle. The total xl ^oi closure phase is 1.33. 



33 - 




Fig. 11. — Similar to Figl9]but showing the triple amplitudes for a Oph. The total for 
triple amplitude is 0.81. 
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Fig. 12. — The xl surface of (3 and inclination for a Oph. The corresponding probabihty 
is high throughout a large range of inclination and P, suggesting high degeneracy between 
the two parameters. The map also indicates the inclination at /3 = 0.25 (i.e., the standard 
model) is well constrained and is nearly equator-on. Since the probability is dominated by 
the degeneracy effects of (3 and inclination, we overplot the xl contours on the map instead 
of confidence intervals. The region enclosed in the green box has Vsini values inside the 
observed range of 210-240 kms~^. The rest of the areas in the map fall outside the observed 
Vsini range and thus can be ruled out, even though they may fit the data better. 




Fig. 13. — Deviation of the apparent Tg// and luminosity from their true values at various 
inclinations. The solid line indicates the standard model (/3 = 0.25) of a Oph. The dashed 
line indicates the /5 = 0.216 model of a Cep. The apparent Tg// and luminosity equal their 
true values at inclination of ~ 54°. a Cep (filled square) is seen very close to this zero- 
difference value, but a Oph (filled dot) is almost at the high end due to its large inclination. 



- 36 - 




9500 9000 8500 8000 7500 7000 6500 

Te« (K) 



AO 



A2 



A5 A7 A9 



50 



40 



^ 30 



20 



10 



1 

: 0.40 Gyr 


0,50 Gyr q gg 


1 1 1 . 

isochrone ' 

evolution track \ 

' -^0^70 Gyr - 
. _ _ _ 0.75 Gyr 






/ _ , Si.m Gyr" 






yd Oph.^ - - - ~ ~ 






"^^r^i - ' ' ' 0:-90 "Gyr- 






_ u - ' ' 1,00 G^i 


2,000Ws 

- 1.9E 

", , , , 1 , , , , 1 , 


920W\s^^^ 

, , , 1 , , , 


■; ':,- ;;-.V:-;-. ;-. , i ,■ 


10000 9500 


9000 8500 

Te„(K) 


8000 7500 



Fig. 14.— Positions of a Cep and a Oph on the HR diagram. The top panel shows a Cep and the stellar models 
with Z~ 0.02. The bottom panel shows a Oph and stellar models with Z ~ 0.014. The dashed red lines indicate the 
evolutionary tracks and the dotted blue lines indicate the isochrones. The filled dots with error bars indicate the true Teff 
and luminosity of the two stars, while the filled triangles indicate their apparent T^ff a-nd luminosity which are dependent on 
their inclinations. The solid lines that go through the poin ts show the positio ns of a Cep and a Oph on the H-R diagram as 
a function of inclination (also called "inclination curves", ICillich et ahlbood) . These curves are more or less p arallel to the 
Zero-Age Main Sequence indicated by the thick solid line at the bottom left of each plot, consistent with those of lGillich et al.l 
liooi). For a 90° inclination, the positions of the stars will be at the lower end of the curve; and for a 0° inclination, the 
stars will be at the higher end of the curves. These plots suggest the inclination of a star can significantly change its apparent 
location on the H-R diagram. 



